Spin-spectral-weight distribution and energy range of the parent compound La 2 Cu04 
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The spectral-weight distribution in recent neutron scattering experiments on the parent compound 
La2CuC>4 (LCO), which are limited in energy range to about 450 meV, is studied in the framework 
of the Hubbard model on the square lattice. We find that the higher-energy weight extends to about 
566 meV and is located at and near the momentum [tt, tt]. Our results confirm that the U/t value 
suitable to LCO is in the range U/t £ (6,8). The continuum weight energy-integrated intensity 
vanishes or is extremely small at momentum [tt, 0]. This behavior of the intensity is consistent with 
that of spin waves, which are damped at [tt, 0]. 



PACS numbers: 78.70.Nx, 74.72.Cj, 71.10.Fd, 71.10.Hf 

It is natural to expect that a greater understanding of 
the physics of the high-T c superconductor undoped par- 
ent compounds, such as La2Cu04 (LCO), will lead to a 
greater understanding of the corresponding superconduc- 
tors. In particular, the less complicated undoped systems 
can provide valuable information on which model Hamil- 
tonians quantitatively describe the cuprates. Improved 
determination of the model Hamiltonians is essential be- 
cause of the many nearby competing phases in the doped 
systems, easily affected by small parameters, which can 
now be seen because of continued improvements in nu- 
merical simulations A decade ago the neutron scat- 
tering experiments on LCO of Coldea, et. al. [|J first 
showed sufficient details of the spin-wave spectrum to 
demonstrate that a simple nearest-neighbor Heisenberg 
model must be supplemented by a number of additional 
terms, including ring exchanges. These terms arise natu- 
rally out of a single band Hubbard model with finite U ft, 
and several detailed studies showed that the spin-wave 
data in the available energy window could be success- 
fully described by the Hubbard model using a somewhat 
smaller value of U/t ~ 6 — 8 than originally thought ap- 
propriate [IHH ■ However, part of the spectral weight was 
deduced to be outside the energy window. 

Recently, improved neutron scattering experiments [f| 
with a much wider energy window of about 450 meV, 
have raised a number of questions. Surprisingly, these 
studies revealed that the high-energy spin waves are 
strongly damped near momentum [tt, 0] and merge into 
a momentum-dependent continuum. These results led 
the authors of Ref. Q to conclude that "the ground 
state of La2Cu04 contains additional correlations not 
captured by the Neel-SWT [spin- wave theory] picture". 
This raises the important question of whether the more 
detailed results can still be described in terms of a simple 
Hubbard model. Here we address this question using a 
combination of a number of theoretical and numerical ap- 



proaches, including, in addition to standard treatments, 
a new spinon approach for the spin excitations 0, 0] and 
density matrix renormalization group (DMRG) calcula- 
tions for Hubbard cylinders [8j4l0|. We show that the 
Hubbard model does describe the new neutron scattering 
results. In particular, at momentum [tt, 0] the continuum 
weight energy-integrated intensity is found to vanish or 
be extremely small. Furthermore, we find that beyond 
450 meV, the spectral weight is mostly located around 
momentum [tt, tt] and extends to about 566 meV, sug- 
gesting directions for future experiments. 



■ Analytical approximation 
t/U expansion 
DMRG, 10x6 cylinder 
DMRG, 10x4 cylinder 



U/t 



FIG. 1: Average single-occupancy: approximate expression 
[1 + tanh(*7/8i)]/2 valid for U/t < 8 (full line), the limiting 
U/t ^> 1 expression [1 — co(8t/U) 2 ] (dashed line), and from 
DMRG numerical results on two different width cylinders. 



The Hubbard model on a square lattice with N S> 1 
sites and periodic boundary conditions reads H = T + 
UD. Here f = -t T,(j,f) E^.t^.o- c v CT + is 
the kinetic-energy operator, D = EjLi ^hj,^^hj,i counts 
the number of electron doubly occupied sites, cj. a and 



c rj,a are electron creation and annihilation operators 
with site index j = 1,.,.,N and spin a =t,i, and 



2 



ct. rrCri.a- The expectation values, 



d = (n rj!t i%. jti ) , (1 - 2d) = ((n rjtt 
m AF = o(( _1 )''(^.t- 



i))n[l-28S]m Q AF ,(l) 



play an important role in our study, following the strong 
evidence that for U > and W-xx) the model ground 
state has antiferromagnetic long-range order . In the 
last expression of Eq. ((TJ) , m AF stands for a mean-held 
sublattice magnetization that does not account for the 
effect of transverse fluctuations and SS accounts for such 
an effect, its value being estimated below. Moreover, 
S» 
S: 



Our study involves the dynamical structure factors, 



S aa (k,w) = 



TV 



E 



3 -ik(rj-iy 



j,i'=i 

dte^*(GS|S?.(t)S^',(0)|GS) , (2) 



where a = a' — x, y, z or a = — and a' = + and be- 
low we consider g — 2. It is straightforward to show 
that the sum rules [l/N]J2 k S aa ' (k) where S QQ '(k) = 
[1/2tt] dio S aa (k, u) involve the average single occu- 
pancy (1 — 2d) and read [(g^s) 2 /4][<5 Q , iQ '+2(5 Q ._J ct ' i+ ](l — 
2d) . In an ideal experiment one sees a transfer of weight 
from the longitudinal to the transverse part such that in- 
dependent of the scattering geometry, the corresponding 
effective spin form factor satisfies the sum rule, 



1 E- 

N 4^ 2tt 



dujS exp (k,uj) = n%2(l - 2d) . 



(3) 



That the coefficient involved is 2(1 — 2d) rather than 
3(1 — 2d) follows from in the experiment one mode being 
perpendicular to the plane and thus silent 

For U/t € (0, 8) the antiferromagnetic long-range order 
may be accounted for by a variational ground state with 
a SDW initial trial state, such as \G) = e~ 9D \SDW) or 



\GB) 



-hT/t p —gD 



state of a simple effective mean-field Hamiltonian, as that 
of Eq. (18) of Ref. 0. For U/t > I such an or- 
der is as well accounted for by a Baeriswyl variational 
state \B) = e~ hT l t \oo) where |oo) is the exact U/t — >• oo 
ground state [Hj]. The coefficients h and g multiplying 
the kinetic-energy and double-occupancy operators, re- 
spectively, are variational parameters. Similarly for the 



trial state \SDW), the relation 4[ 



"AFl 



= (l-4d) holds 



for \G), \GB), and \B), while the function d = d(U/t) is 
state dependent. This gives d = j[l — 4[ttt,5i^] 2 ], consis- 
tent with d not being affected by transverse fluctuations. 

The evaluation of the ground-state energy for \GB) 
and \B) is for N 3> I an involved problem. Here we resort 




(3/4,1/4) 



FIG. 2: Upper panel: Spin-wave excitation spectrum along 
Brillouin-zone special directions as specified in Ref. @] . Lower 
panel: Spin-wave intensity as obtained from the poles of the 
susceptybility (see text). Experimental points from Ref. [f|. 



to an approximation, which corresponds to the simplest 
expression of the general form E/N = T^qu + Ud where 
To = — ^4 t compatible with three requirements: the rela- 
tion d — |[1 — 4[m^ F ] 2 ]; the occurrence of antiferromag- 
netic long-range order for the whole U/t > range; and 
the lack of a linear kinetic-energy term in U for U/t <C I . 
Brinkman and Rice found qjj = 8d(l — 2d) for the original 
Gutzwiller approximation [13| , which is lattice insensitive 
and thus does not account for the square-lattice antiferro- 
magnetic long-range order. The simplest modified form 
of the quantity qu such that the three conditions are met 



a\ +) d 



(l-2ri) 



- a 2 



A3 



where 4[m 

:±) 



\SDW), where | SDW) is the ground a 2 



is qu — V8t7 "l 

behaves as 4[m AF } 2 = U/8t for U/t < 1, U| 

(1 — [tt 2 /2a^}), and as — aa(U/8t) a function given 



AFl 

2 ±4, 



by a 3 = [a{^/8}{l - tanh([l//8t][(4 + a^')/a^'])} for 
U/t G (0, 8) and a 3 = -c [ir/2} 2 [8t/U] for U/t > 1. Here 
Co = [a/4 + 1/8]/ 2 = 0.1462 and the corresponding es- 
timate a = 0.6696 is that of the Heisenberg-model stud- 
ies of Ref. [HI]. Minimization of the obtained ground- 
state energy with respect to d leads indeed to d = 
j[l — 4[to^ f ] 2 ]. For U/t <C 1 such a qu expression yields 
to second order in U/t, E/N = T + Ud- [1 /8n 2 ][U 2 /t] 
where [1/8tt 2 ] rj 0.0127. The coefficient -0.0127 
is that also obtained by second-order perturbation the- 
ory [3] • For U/t ^ 1 one recovers the known result 
E/N = — Aco[8t 2 /U] (ll| . so that our approximation 
agrees with the known limiting behaviors. 



(+h /J-h 
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For U/t G (0,8), we find that 4[rn^] 2 w tanh(t//8t) 
gives quantitative agreement for the (1 — 2d) dependence 
on U/t with both our numerical DMRG calculations (see 
Fig. 1) and the numerical results for the states \G) and 
\GB) (see Fig. 4 of Ref. d). In the DMRG calcula- 
tions, two different circumference cylinders were simu- 
lated as a function U /t, with open boundary conditions 
in x and periodic in y, and the double occupancy mea- 
sured in one of the middle columns. A maximum of 
rn = 4000 states were kept, with an accuracy of ~ 10~ 4 
in (1 — 2d) for the 10 x 4 system for the least accurate 
smaller U /t values, and about 10 -3 for the 10 x 6 system. 
We find that the value of (1 — 2d) is relatively insensitive 
to cluster size, and these cluster sizes are representative 
of 2D behavior 15]. For [J/i > 1 we find the behavior 
4[m° F ] 2 » e -2co (8t/C/) 2 for the gtate SQ that _ 

2d) = [l+tauh(f7/8i)]/2 for approximately U/t < 8 and 
(l-2d) = [l-c (8t/U) 2 ] for U/t > 1. Furtherm ore, the 
states | SD W) and \B) give m AF « m AF = |\/l — Ad, 
with an improved d = d(U/t) dependence for the latter, 
whose m AF magnitudes are provided below in Table Q] 
for several U/t values. The states \GB) and \B) have 
m A p and m AF sublattice magnetization numerical val- 
ues very close to those given by Eq. (JTJ with SS « d and 
55 « d + | [1 - m HAF /m° HAF } 
Here fn {) HAF 



d + 0.197, respectively. 
1/2 and ttihaf ~ 0.303 is the Heisenberg- 
model sublattice magnetization magnitude [ill ]. Some 
to J|J? magnitudes are given below in Table HI along with 
those of m l AF er = (1 — 2d)[mHAF/m ( j IAF ] to^f, which 



'AFi 



we define for U/t > and for U/t 3> 1 becomes m AF . 
Probably m AF er is closest to the exact tuaf, while m§J? 
is consistent with our use of the RPA. 



To study the coherent spin-wave weight distribu- 
tion and spectrum, we have calculated by RPA 
the transverse dynamical susceptibility x (kj r ) 



ik -( r i- r i)(5-(r J -,r)S+(r J v,0)), where r 



denotes the imaginary time in Matsubara formalism (we 
shall take the zero temperature limit). In the case of the 
present model, the coherent spin-wave spectral weight in 
units of fi% is Z d 2(l - 2d). Here Z d = I - [2/JV(l - 

2d)] SkS„'^i, K^'Pk l^^)! 2 where the sum over energy 
eigenstates excludes those that generate the coherent 
spin-wave weight, \v) = \v,u)(k)). In the [//i — > oo limit, 
Zd may be identified with the corresponding Z d — Z c Z x 
factor of the Heisenberg model on the square lattice. 
According to the results of Ref. [l(|, Z c w 1.18 and 
Z x w 0.48, respectively, so that Z d w 0.57. The lim- 
iting values Z<£ = 1 for U/t -> and Z d ss 0.57 
for U/t oo and the approximate intermediate value 
Zrf w 0.65 at U/t — 8 4] are recovered as solutions of 

the equation Z d = e Zd tanh ("V / S) ) wn i cn j s usec j here 
for finite {/ /t. In the thermodynamic limit the upper- 
Hubbard band processes generate nearly no spin weight. 
Hence the longitudinal (elastic contribution to) spectral 
weight within fi%2(l — 2d), Eq. ([3J, is in units of n 2 B 




FIG. 3: The energy-momentum space limits of the spin S — 1 
excited states spectrum for U/t = 6.1, t = 295 meV, and k x 
and k y in units of 2n. States whose energy is for a given k 
lower than that of the intermediate spin-wave sheet as well as 
those of any energy and equivalent momenta [0, 0] — [0, 2-k] = 
[2-7T, 0] = [2-7T, 2n] do not contribute to the spin spectral weight. 



approximately given by ~ A^uiaf) 2 and the spin- wave 
intensity reads W sw = Z d [2(1 - 2d) - 4(m A F) 2 }- The 
GA+RPA method used in Ref. [4[ accounts for the 
quantum fluctuations that control the longitudinal and 
transverse relative weights, so that the spin-wave inten- 
sity factor is Zd- The RPA used here leads to a similar 
spin-wave intensity momentum distribution but its ex- 
_ perimentally determined factor Z e d xp < Z d is such that 



Z e d xp 2(l-2d). 



W SW = Z d [2(1 - 2d) - 4(m AF ) 2 ] 

Within the two-sublattice description the susceptibil- 
ity becomes a 2 x 2 tensor and the above original suscep- 
tibility x h (k, t) is the average of its four elements. Af- 
ter Fourier transforming to (k, iu) space, within the RPA 
the susceptibility tensor obeys x RPA = [l — Ux^^ 1 ^ * 1 - 
Here x'' '' denotes the susceptibility tensor in the nonin- 
teracting system. x RPA has a pole ioj = oj(k) obtained 
from the equation Dot [1 — Ux^] = 0, which provides 
the dispersion relation w(k) for the spin waves. It has 
been shown in Ref. Q that an excellent agreement with 
the spin-wave spectrum from Ref. 0] is achieved. In Fig. 
[2] upper panel we show a fit to the more recent experi- 
mental data of Ref. (full line) alo ng with the results 
from the si fermion method reported below (dotted line) 
for U/t = 6.1 and t = 295 meV. Provided that the t 
magnitude is slightly increased for increasing values of 
U/t, agreement with the LCO spin- weight spectrum and 
distribution can be obtained for U/t G (6,8). For U/t 
values smaller than 6 (and larger than 8), the spin- wave 
dispersion between [ir/2, 0] and [tt/2, 7r/2] has a too large 



4 



energy bandwidth (and is too flat) for any reasonable 
value of t. From the residue of the spin- wave pole the 
susceptibility coherent part reads, 



l=±l 



Res [x~ + (k,/o;(k)] 
iu> — I uj(k) 



(4) 



The measured intensity is Jsw(k) = Tr[S xx (]i) + 
S m (k)} = 7r5-+(k) [17]. In Fig. flower panel, we plot 
the corresponding RPA spin- wave intensity, Isw(k) = 
— [7r/2]Z^ p Res [x h (k, w(k)]. The good agreement with 
the experimental data, specially near the point M, re- 
produces the theoretical results of Ref. [6]. It is here 
obtained for the value Z^ xp ps 0.49, which corresponds 
to the choice tuaf = m AF = 171 af(^ ~ %d) such that 
SS ps d. As in Ref. Q , it shows disagreement around the 
X point, which here probably stems from the RPA. 

To derive the spectrum of the spin 5=1 excited states 
producing the inelastic form-factor spectral weight, we 
use the spinon operator representation of Ref. 15]. The 
ground-state spin degrees of freedom are described by 
a full si fermion band, which coincides with an anti- 
ferromagnetic reduced Brillouin zone. Each si fermion 
is a spin-singlet two-spinon composite object. The spin 
5 = 1 excited states of momentum k = [it, tt] — q — q' 
are generated upon breaking one si fermion spinon pair, 
which leads to the emergence of two holes of momenta 
q and q' in the otherwise full si band. The inelastic 
coherent spin-wave spectrum is generated by processes 
where q points in the nodal direction and q' belongs 
to the boundary of the si band reduced zone [§]. The 
remaining choices of q and q' either generate the in- 
elastic incoherent continuum spectral weight or vanish- 
ing weight, respectively. The studies of Ref. [B[ are lim- 
ited to the spin-wave spectrum. Here we consider the 
energy-momentum space domain of all above spin 5=1 
excited states, which is represented in Fig. i3l (A sim- 
ilar spectrum is obtained for the values U/t — 8.0 and 
t = 335 mcV of Ref. ,4j].) The intermediate sheet refers 
to the spin-wave spectrum. For each k, states of energy 
lower than the latter spectrum do not contribute to the 
form-factor weigh. Furthermore and consistent with the 
first-moment sum rules of an isotropic antiferromagnet, 
no and nearly no weight is generated by states of any 
energy and momentum [0, 0] and near [0, 0], respectively. 
Unfortunately, the method of Ref. Q does not provide 
the detailed continuum weight distribution. However, it 
is expected that, similarly to the Heisenberg model case 
0] , its energy- integrated intensity follows the same trend 
as the spin- wave intensity. Analysis of Fig. [3]reveals that 
for momentum [tt, 0] there are no excited states of energy 
higher than the spin waves. Thus at momentum [tt, 0], 
the continuum weight distribution energy-integrated in- 
tensity vanishes or is extremely small, due to si band 
four-hole processes. Given the expected common trend 
of both intensities, this is consistent with a damping of 



U/t 


6.1 


6.5 


8.0 


10.0 


W T 


1.643 


1.671 


1.762 (1.778 [4]) 


1.848 (1.846 [4]) 


Wsw 


0.808 


0.799 


0.761 


0.714 


W<450 


1.0/1 


i.oyo 


l.DDo 


1. / OU 


W>450 


0.072 


0.078 


0.099 


0.118 


m AF 


0.401 


0.410 


0.436 (0.43 [4, .12]) 


0.461 (0.456 [4J) 




0.329 


0.342 


0.384 (0.39 [12]) 


0.426 




0.200 


0.207 


0.233 


0.258 



TABLE I: Several spectral weights in units of /j% as defined in 
the text and the sublattice magnetizations as calculated here 
for several U/t values and some results from Refs. [H1I3]- 



the spin- wave intensity at momentum [tt, 0], as observed 
[6j but not captured by the Fig. [2] RPA intensity. 

In Table U we provide the sublattice magnetization 
magnitudes along with our calculations of the follow- 
ing integrated spectral weights (in units of fi%): the 
total weight, Wt = 2(1 — 2d); the spin-wave weight, 
Wsw = Zf p 2{l - 2d)- the weight for energy Hlu <450 
meV, VF<450 = Wsw/0.71 + 4(m^f ) 2 ; and the weight 
W>450 = [Wt — lF<45o] for energy Huj > 450 meV. W^o 
is derived by combining our theoretical expressions with 
the observations of Ref. @ that for the energy range up 
to about 450 meV, 71% and 29% of the weight corre- 
sponding to the inelastic response comes from the coher- 
ent spin- wave weight and incoherent continuum weight, 
respectively. Our prediction for VF<45o varies between 
1.6(4% for U/t ps 6.1 and 1.7/4% for U/t ps 8.0, in agree- 
ment with the experimental value 1.9 ± 0.3 /x% reported 
in Ref. 0. From our above analysis, the small weight 
1^>450 ~ 0.1 (i% is expected to extend to about 566 meV 
(see Fig. [3]), mostly at and around the momentum [n, tt]. 

Our result that at momentum [tt, 0] the continuum 
weight energy-integrated intensity vanishes or is very 
small is consistent with a corresponding damping of the 
spin- wave intensity at [tt, 0], as observed in the recent 
experiments of Ref. @. We suggest that future LCO 
neutron scattering experiments scan the energies between 
450 meV and 566 meV and momenta around [ir, tt]. 
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